Nucleon-Meson Coupling Constants and Form Factors 
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We demonstrate the calculation of the coupling constants and form factors required by effective 
hadron lagrangians using the quark model. These relations follow from equating expressions for 
strong transition amplitudes in the two approaches. As examples we derive the NNm nucleon- 
meson coupling constants and form factors for m = ■n,r],ri' ,cy,ao,uj and p, using harmonic oscillator 
quark model meson and baryon wavefunctions and the ^Po decay model; this is a first step towards 
deriving a quark-based model of the NN force at all separations. This technique should be useful in 
the application of effective lagrangians to processes in which the lack of data precludes the direct 
determination of coupling constants and form factors from experiment. 

PACS numbers; 12.39.Jh, 13.75.Gx, 21.30.Fe, 24.85. +p 



I. INTRODUCTION 

Effective hadronic lagrangians are widely used the de- 
scription of the interactions of hadrons. In this method 
a distinct quantum field is introduced for each relevant 
hadronic species, and interactions are assumed for these 
fields that are consistent with known symmetries and 
conservation laws. Although a more physically justified 
description of hadron scattering amplitudes would em- 
ploy momentum-dependent coupling constants (form fac- 
tors), it often suffices near threshold to assume pointlike 
"hard" hadronic vertices with fixed coupling constants. 
This approximation may be relaxed by incorporating 
hadronic form factors as a power series of gradient inter- 
actions. Unfortunately, this typically leads to nonrcnor- 
malizable ultraviolet divergences, which requires the de- 
termination of many coupling constants from experiment. 

When using effective lagrangians one typically ignores 
the existence of hadronic substructure (quarks and glu- 
ons), and determines each coupling constant in the ef- 
fective lagrangian from experiment. Although this pro- 
cedure is feasible in experimentally highly constrained 
problems such as the NN interaction, in other processes 
with little data the coupling constants are simply as- 
signed plausible values. In such cases it would be very 
useful to have numerical estimates of the coupling con- 
stants. Since these coupling constants and form factors 



are determined by the underlying QCD degrees of free- 
dom, one may evaluate them directly in terms of the 
interactions of quarks and gluons and the substructure 
of the hadronic bound states. 

In this letter we will investigate this relation between 
effective lagrangian couplings and quark model bound 
state hadron wavefunctions. The specific cases we con- 
sider are the meson-nucleon couplings, which are usu- 
ally fitted to data in meson exchange models of the NN 
force. These NNm vertices are chosen as our initial exam- 
ples in part because they are among the most important 
hadronic couplings for nuclear physics applications, and 
also because the NNtt coupling constant ^nntt ~ 13.5 is 
the best determined strong coupling in hadron physics. 
It is clearly of interest to determine whether this exper- 
imental coupling constant is consistent with the value 
predicted using quark model hadron wavefunctions. 

We regard this computation as the first step in an at- 
tempt to compute the strong interactions of nucleons at 
all separations from microscopic quark-based models. In 
this approach the short ranged NN interactions arc dom- 
inated by the quark-gluon forces encountered by overlap- 
ping nucleon quark wavefunctions, and the long ranged 
NN interactions are mediated by t-channel meson ex- 
change. The meson-nucleon couplings in t-channel me- 
son exchange, which we compute here, are themselves 
determined by the overlap of meson and nucleon quark 
wavefunctions. 
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II. TECHNIQUE 

Our method for defining eff'ective hadronic coupling 
constants and form factors is to equate specific hadron 
emission ampfitudes predicted by the efi^ective hadron la- 
grangian to the corresponding decay amphtude in the 
quark model. Applied to NNm couplings in particular, 
we require that 

{Nm\H,ff\N) = {N{q'')m{qq)\Hdecay\N{q'')) . (1) 

This general approach was used previously by the Orsay 
group of LeYaouanc et al. 0, who considered pmv, NNtt 
and NNp couplings; the relation between this earlier ref- 
erence and our results will be discussed. A similar quark 
model approach was more recently applied to the deter- 
mination of HQET electroweak form factors by Isgur et 
al. 

Ideally one would impose this relation using a rela- 
tivistic quark model, in which case there would be no 
difficulty in identifying the effective lagrangian matrix 
element with the quark model result. As the nonrel- 
ativistic quark model formalism is much better estab- 
lished, we will instead use nonrelativistic quark model 
wavefunctions, and apply our defining relation Eq.(^ be- 
tween matrix elements near threshold. Since the quark 
model wavefunctions we use are nonrelativistic, the form 
factors we extract by equating quark model and (rela- 
tivistic) field theory matrix elements have an ambiguity 
in how we relate results in different frames. In this ini- 
tial study we will simply assume a particular frame, the 
"decay frame" of a rest initial hadron, as in the earlier 
work of the Orsay group The complication of im- 
plicit frame dependence will be investigated in detail a 
subsequent study. 

It is important to note that these frame ambiguities are 
resolved in the weak binding and nonrelativistic limits, in 
which the quark model states approach the transforma- 
tion properties specified by the Lorentz group. Thus our 
quark model predictions of form factors are best justified 
at zero three-momentum, where we define the coupling 
constants. Of course the model assumption is that the 
predicted form factors remain useful estimates at moder- 
ate nonzero momentum. 

The specific quark model strong decay Hamiltonian we 
use is the model which has seen very wide 

application, and is known to give reasonably accurate 
numerical results for both meson 0, IS ^^'^ baryon Q 
strong decays. 



A. NNtt, NNr; and NNrj' 

As a first example we consider the NNtt coupling. We 
assume the standard pseudoscalar effective lagrangian 

-Cnntt = -ig-NN^'^ToT'^ ■ cj)-^ , (2) 
where 5* = [ipp,ipn] is an isodoublet of Dirac nucleon 
fields, the r are isospin Pauli matrices, and (^^r is an 



isotriplet of pion fields, [07r+i '/'ttC'i '/'tt-]- Specializing to 
the transition j3(+l/2) p(-|-l/2) tt" (an initial spin-up 
proton at rest going to a spin-up proton with momentum 
— P and a recoiling tt*^ with momentum -l-P), we find the 
matrix clement 



{pAHeff\p)/5{^S-^^)^h^ 



Mr. 



«5nn7t 



E-p ^(27r)32K ^^2Mp{Ep + Mj 



P cos(6l) 



(3) 



We similarly evaluate the quark model matrix element 
of the ^Po decay model Hamiltonian for p{+l/2) 
p{+l/2) TT^, using the techniques and Gaussian baryon 
and meson wavefunctions given in Refs.0,|3- The result 
is 



{piq^)n°iqq)\Hdecay\piq^))/SiPf 



10 (l + rV4)r3/2 Pcos(6i) 
a3/2 



h 



97r3/4 (1 + ^2/3)5/2 



exp 



1 (l + 5rVl2) 
6 (l + r2/3) V" 



(4) 



where a and f3 — a/r are the baryon and meson Gaus- 
sian quark wavefunction inverse length scales, and 7 is 
the dimensionless '^Pq model qq pair production am- 
plitude. This result is consistent with the earlier re- 
sult of LeYaouanc et al. 0, given the parameter re- 
lations between that reference and the current work, 
(i?N = 1/a), {R^ = 1//3) and {-/LeY = 724^7)- Typi- 
cal values for these parameters found in the literature are 
a = 0.25-0.4 GeV, (3 = 0.3-0.4 GeV, and 7 = 0.4-0.5. 
Here we fix 7 = 0.4 (taken from our extensive studies of 
light and charmed meson strong decays 1, "F, ) , and 
show numerical results for this range of a and /3. 

The NNtt coupling constant wc find by equating these 
expressions at threshold is 

24.5 3/4 (l + rV4)r3/2 Mpm'J^ , , 
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(l + r2/3)^/2 a3/2 



(We have suppressed an overall phase factor of -l-i, which 
is normalization convention dependent.) The numerical 
value of this coupling constant as a function of a and /3 is 
shown in Fig^ Evidently the range of typical wavefunc- 
tion length scale parameters a and /3 leads to a factor 
of two variation in the theoretical (/nntt! it ranges be- 
tween 7.0 — 12.2 for the parameters shown in the figure. 
The experimental value of ~ 13.5 evidently corresponds 
to values of a and (3 near the lower end of their respec- 
tive ranges, provided that we fix the qq pair production 
amplitude at the meson decay value of 7 = 0.4. 

Since the "experimental" value of ^nntt is not actually 
based on a direct observation of pion emission, it is pru- 
dent to carry out an independent calculation of a closely 
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FIG. 1: The theoretical pion-nucleon coupling constant gNNyr 
(Eq|3 as a function of the quark model Gaussian wavefunc- 
tion length scales a (baryons) and /3 (mesons). The range of 
values of a and /3 typically found in the quark model litera- 
ture is indicated (see text). The ^Po model qq pair production 
amplitude 7 = 0.4 was taken from studies of meson decays. 
The experimental pnntt ~ 13.5 is also shown. 



related decay process that does involve a detected pion 
in the final state. The transitions A Ntt are useful 
in this regard because the matrix elements are related to 
the NNtt coupling by SU(6) symmetry, assuming identi- 
cal spatial wavcf unctions. Specializing to P7r+, 
the quark model result for the total width is 



r(A^ 



Ma 



P 

- ) exp 

a 



2 1/2^ 3(1 
7 33^ (1 



5rVl2) 



-rV4)2 



1(1 

3 (l + r2/3) 



(6) 



This theoretical quark model decay rate is shown in Fig|3 
for the same range of wavefunction parameters a, f3 as 
^nntt in FigE and the experimental width of 110 MeV 
is also shown. Evidently the parameter constraints due 
to (jfNNir and r(A Ntt) are approximately consistent. 

Alternatively, since the dimensionless pair production 
amplitude 7 represents poorly understood nonperturba- 
tive physics, it may have a different strength in baryon 
decays than in meson decays. In Fig|31we show the ra- 
tio of theory to experiment for r(A'''+ — > p7r^) and ^nntt 
for our "standard" quark model baryon and meson wave- 
function length scales a = P = 0.4 GeV, varying the pair 
production amplitude 7. Evidently with these wavcfunc- 
tions a value of 7 « 0.7 is favored for pion emission from 
light baryons, whereas 7 ~ 0.4 — 0.5 gives the best de- 
scription of light meson decays. 

Nucleon couplings to the other ground state pseu- 
doscalars {rj and 77') are simply related to g-MN-Tr, provided 
that we assume identical spatial wavefunctions and pure 
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FIG. 2: The theoretical ^Po quark model r(A++ p7r+) 
width (Eqinj as a function of the wavefunction length scales 
a and /3, for a pair production amplitude of 7 = 0.4. The 
experimental width is also shown. 




FIG. 3: The ratio of theory and experiment for r(A+^ ~* 
piT^) and (?NN7r given our standard quark model wavefunction 
length scales a = f3 = 0.4 GeV, but with the pair production 
amplitude 7 variable. 



qq states. Given the effective lagrangian 

-CnNjjC) = -*5'NN),(')^'>'5* '?^'7(')' 



(7) 



and taking the \ri) and \ri') flavor states to be the max- 
imally mixed linear combinations (\nh ) ± \ss))/V2, the 
NN77(') coupling constants and form factors are related 
to the NNtt results by 

3 /"^J7(')\^/^ 
<?NN,(') = ^T7J:^(^^j (8) 

If we use .gNNTT = 13.5 as input, this gives (/nn?? — 11-5 
and (7nnjj' = 15.3. Although these appear to be rather 
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large NN7y(') couplings, their effect on NN scattering is 
suppressed by the larger 77 and rj' masses in propagators 
and in the Xj \J1E external line normalizations. 



B. NNcr 

The NNcr coupling may be the most important 
nucleon-meson coupling in all of nuclear physics. In 
meson exchange models the exchange of a light scalar 
1=0 "sigma" meson is held to be the dominant mech- 
anism underlying the intermediate ranged attraction, 
which is responsible for nuclear binding. The balance 
between this attraction and the short distance repulsion 
in the nuclear core determines the equilibrium density 
of no ~ 0.16 nucleons/fm^ in bulk nuclear matter. (Al- 
though pions are lighter and hence longer-ranged, and 
the (7NN7r coupling constant is quite large, the fact that 
pions are emitted in a relative P-wave suppresses their 
contribution to the near-threshold interactions of nucle- 
ons.) 

Although the intermediate ranged attraction plays a 
crucial role in nuclear physics, the a meson itself is a 
dubious concept in meson spectroscopy. In 1=0 tttt S- 
wave scattering one sees a very broad positive (attractive) 
phase shift, which if interpreted as an s-channel reso- 
nance would imply a mass of ca. 1 GeV and a comparably 
large width. There are arguments from the quark model 
against a qq state with these properties; for example, an 
1=0 O^"*" nn resonance (n = u, d) at this mass would 
have a rather large two-photon width of F-y^ k, 2 keV, 
and no such state is evident in 77 — > tt^tt*'. There is 
instead evidence in this reaction for a moderately broad 
scalar enhancement near 1.3 GeV, with about the ex- 
pected two-photon width [^; this broad /o(~ 1300) is 
often identified with the 1=0 0++ nn quark model state. 

The explicit cr meson included in meson exchange mod- 
els has been explained as a parametrization of "correlated 
two-pion exchange" , so that its fitted strong coupling to 
NN and low mass need not correspond to the properties 
of a physical P-wave nn meson. The picture of the "sigma 
meson" as a strongly mixed {nn)-Tn: state is supported 
by the large coupling predicted between these channels 
in the ^Pq model; the analogous S-wave kaon system is 
discussed in Ref . |11| . 

We can test the plausibility of sigma meson exchange 
models of the intermediate ranged NN attraction by cal- 
culating the NNcr coupling for a pure nn a meson, us- 
ing the same techniques we applied above to the NNtt 
system. If the sigma is dominantly a physical nn scalar 
meson, we would expect approximate agreement between 
the calculated and fitted (^nnct coupling constants. If the 
sigma is instead a parametrization of two-pion exchange, 
agreement between the theoretical and fitted coupling 
constants would be fortuitous. 

The calculation of the NNcr coupling differs from the 
NNtt case through the meson spin, space and isospin 
wavefunctions and the effective lagrangian. We assume 



the form 



-NNo 



-5NNcr^'*< 



(9) 



In our quark model description the |cr) is the usual 1=0 
\nn) flavor state {\uu) + \dd))/\/2 times the |J,L, S) = 
|0, 1, 1) angular momentum state 

-i=(|l,+l)|l,-l) - |1,0)|1,0) + |1,-1)|1,+1)) 

where the basis states are |L, Lz)|S, Sz). The P-wave 
momentum space qq wavefunctions are given in Ref.Q. 
On equating the effective lagrangian and quark model 

hfi — {pa\H\p) matrix elements, we find 



5NNt 



= 7 25/231/2^3/4 



r.5/2 



1/2 



(l + r2/3)5/2 ai/2 



(10) 



(A normalization convention dependent overall phase of 
(—1) in our result is suppressed here.) The NNcr form 
factor is the quadratic {l + [{l+r'^ / 4) /^l+r"^ af) 
times the Gaussian found in the NNtt case in Eq.l^J. 

The numerical (/nno- predicted by Eg. 1)10(1 is shown in 
Fig^as a function of a and /3 (assuming nia — 500 MeV). 
Evidently a value in the range 3 — 7 is predicted by the 
quark model given this m^, with (Jnnct ~ 5 preferred. 




FIG. 4: The theoretical nucleon - sigma meson coupling con- 
stant gNNo- fEallOfl. 

Although it is of great interest to compare our calcu- 
lated (jNNo- coupling constant with the values reported in 
meson exchange model fits to NN scattering data, there 
is unfortunately no single consistent value reported for 
(7NNcr in these models. The three best known meson- 
exchange models of the NN force in the literature are 
the Paris [H [H, Nijmegen [H [H [3 and Bonn 
[TtI models, and their NNm couplings are given in Ta- 
ble HJ together with our quark model results. The Paris 
model did not consider a a meson. In the recent "CD- 
Bonn" model [1^ , different (7nN(t coupling constants and 
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a masses are assumed in different NN channels; in the 
1=0 ^Si NN channel a a mass and coupling constant of 
rria = 350 MeV and ^nno- ~ 2.5 are used, whereas in 1=0 
^So, ma = 452 MeV and (7nn<t ~ 7.3 are used. The ^nno- 
coupling is allowed to vary with L and I in the L > NN 
channels, and ranges from 1.9 to 9.9 (see Table XVI of 
Rcf. In the Nijmegen model |0| a larger value of 

5NNcr = 17.9 is quoted. Thus, although our quark model 
result (7NNcr ~ 5 is similar to the mean S-wave NN value 
quoted in the CD-Bonn model, the scatter in the fitted 
values of this parameter precludes an accurate compari- 
son between theory and experiment at present. 

One may similarly evaluate the quark model prediction 
for the NN coupling of the aq 1=1 partner of the a. Given 
the effective lagrangian 



L 



we find 



NNao 



5NNao — o \ 

3 V nin 



-gNNao*^' * • 



1/2 



.9nno 



(11) 



(12) 



Although 1=1 scalar exchange contributes a somewhat 
smaller amplitude to NN scattering than 1=0 exchange, 
it may nonetheless be possible to test for the presence of 
both of these scalar meson exchange amplitudes through 
their interference, for example by comparing the 1=0 and 
1=1 S-wave NN scattering amplitudes discussed above. 



C. NNtj and NNp 

The NNV couplings are interesting in that the short- 
ranged repulsive core in the NN interaction has previ- 
ously been attributed to vector meson exchange, and the 
existence of the lo meson was regarded as support for 
this picture. This mechanism now appears less plau- 
sible, since the short range of vector meson exchange 
{R ~ 1/mij ~ 0.25 fm) implies extensive overlap of the 
NN quark wavefunctions. 

Evaluation of the NNV couplings and form factors uses 
the same procedure as the scalar and pseudoscalar cou- 
plings discussed above, although there are complications 
due to the presence of two form factors and the non- 
transverse components of the vector field. 

As above we assume a term in the effective lagrangian 
for each coupling. For NNw this lagrangian is 



L 



4Mn 



^cT^.u^'P^.u] (13) 



where F^i^ = — d^cu^. We then equate near- 

threshold Hamiltonian matrix elements hfi found from 
this effective lagrangian to the corresponding "^Pq decay 
model matrix elements. There is a complication in re- 
lating the relativistic effective lagrangian and nonrela- 
tivistic quark model matrix elements; we find that one 
must assume a vector meson polarization vector and four- 



equate these expressions. The NNw 7^ and form fac- 
tors may be separated by equating hfi matrix elements 
with different spin states. The transitions p(-|-l/2) — > 
p(+l/2) cj(0) and p(+l/2) ^ p(+l/2) are useful 

in this regard, since they receive contributions from only 
the 7^ and cr^j/ terms respectively. The resulting form 
factors are proportional to the NNtt result Eq.©, since 
they involve the same Gaussian overlap integrals. The 
NNw coupling constants (and form factors) satisfy the 
relations 



9 /m^\ 1/2 
5 V m^r / 



SNNtt 



and 



= -3/2 



(14) 



(15) 



The NN/9 form factors, defined through the generalization 
of the NNo; effective lagrangian to an 1=1 p meson, 



-Cnnp = -5NNp (^^'7M^*-/0M-^^*f^M^^*-FM^j (16) 

(where — d^p^ — dyPf^) are related to the NNw results 
by 



1 /m„\i/2 



and 



1 / mp \ 
3 V m„, / 



Kp = -f 3/2 . 



(17) 



(18) 



momentum of the form e,, 



(0, e) and = (0, q) to 



The NNp vector (7^) coupling constant was previously 
evaluated by LeYaouanc et al. fl'|. Their Eq.(3.13) for 
.9NNp is consistent with our Eas. (|5ll7|l . provided that i) 
their factor of SR'^R'^ is actually 3R^Rp (their result 

as written is dimensionally incorrect), ii) their factor of 

3/2 1/2 
TOp should instead be nip ttin, analogous to their (?NNir 

result, and iii) the factor of in their p effective la- 
grangian Eq.(2.17) should be r, as in their tt effective 
lagrangian Eq. (2. 12). LeYaouanc et al. did not evaluate 
the (T^y term, and did not consider the NNw case. 

It is interesting to compare our theoretical NNV cou- 
plings with the fitted values found in meson exchange 
models of NN scattering. If the '^Po model is reason- 
ably accurate in describing the coupling between non- 
strange baryons and vector mesons (which is currently 
being tested at TJNAF in their search for missing baryon 
resonances decaying to No; and Np), and the meson ex- 
change models are correct in assuming that vector me- 
son exchange is the dominant mechanism underlying the 
short-ranged NN force, wc would expect to find approx- 
imate agreement between these couplings. 

The fitted NNV couplings found in the three well- 
known meson exchange models are given in Table H to- 
gether with our quark model results. Evidently we do not 
find good agreement. Note in particular that the fitted 
strength of the dominant NNlu 7^ coupling in the meson 
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exchange models is about a factor of 2 smaller than the 
quark model prediction. The ratio of the NNp to NNw 7^ 
couplings is rather similar in the two approaches; the me- 
son exchange models quote a ratio of « 0.2 — 0.4, whereas 
the theoretical ratio (an SU(6) symmetry factor rather 
than a detailed test of the quark model predictions) is 
+1/3. Although the ratio "k" of magnetic (<t^i/) to vector 
ilfj.) couplings does not yet appear to be well determined 
for both vectors in the meson exchange model fits, there 
does appear to be agreement that \Kp\ » \ku\. This 
is inconsistent with our quark model prediction of equal 
magnitudes for these NNV "strong magnetic" couplings, 
K(uj.p) = (— ,+)3/2. Since these ratios are simple SU(6) 
factors and do not involve uncertainties in the spatial 
wavefunctions, this disagreement may imply that vector 
meson exchange is not the dominant short ranged NN in- 
teraction mechanism. This will be addressed in detail in 
a future study of the NN scattering amplitudes and phase 
shifts due to meson exchange, augmented by quark model 
constraints on the nucleon-meson vertices. 



III. SUMMARY AND CONCLUSIONS 

In this paper we have developed a formalism for deter- 
mining hadron strong vertices and form factors, "three- 
point functions" , in the context of the quark model. We 
apply this approach to the evaluation of meson-nucleon 
vertices, several of which are important in meson ex- 
change models of nuclear forces. The quark model ex- 
pression we find for the NNtt coupling confirms an earlier 
Orsay result. With standard quark model parameters, 
this (^nntt about half the experimental value. Our quark 
model expression for the theoretical gNNo- strong coupling 



of nucleons to scalar mesons is a new result, and is numer- 
ically similar to the isospin-mean fitted NN S-wave value 
in the CD-Bonn model. Our quark model result for the 
NNp vector (7^) coupling is consistent with an earlier 
Orsay result (after correcting typographical errors), al- 
though we find a nonzero magnetic (cr^;y) coupling. The 
strengths of the fitted NNV 7^ couplings in meson ex- 
change models are rather smaller than our quark model 
predictions. The NNV a^^ couphngs are also not in good 
agreement with quark model predictions, although they 
may not be well determined in the current fits to NN 
scattering data. 

In future we plan to carry out calculations of the NN 
phase shifts predicted by meson exchange models, assum- 
ing quark model constraints on the NNm couplings and 
form factors as derived here. This should allow a deter- 
mination of the sensitivity of the data to parameters such 
as the gNNoj/SNNp and Kuj/Kp ratios, for which we have 
definite quark model predictions. 
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TABLE I: A summary of NNm coupling constants. Our calculated values are shown in the middle columns, and fitted 
assumed values in the meson exchange model literature are shown at right. Values in square brackets were fixed input. 



Coupling 


This ref." 


This ref." 


Paris [12, 13] 


Nijmegen [15] 


Bonn [17] 


QNNn 


7.1 


[13.5] 


[13.3] 


13.3 


[13.1] 


ffNNr; 


6.0 


11.5 




9.8 






7.9 


15.3 


10.5 


<7nN(t(~500) 


5.0 




17.9 


(2.5; 7.3)" 


<?NNao(~1300) 


2.7 


3.3 


QNNu, (7m) 


30.2 


12.2 


12.5 


15.9 


gNNp/gNNoj (7m) 


+1/3" 


0.43" 


0.22 


0.20 




-3/2 


-0.12 


0.66 







+3/2 




6.6 


6.1 



"Assumes "standard" quark model parameters a = (3 = 0.4 GeV 
and 7 = 0.4 (see text). 
''Assumes ^nntt = 13.5. 

'^This "CD-Bonn" model introduces different flNNo- coupling con- 
stants for (1=0; 1=1) NN channels, which would not be expected 
for an isosinglet a. In addition these gNNo- couplings and the a 
mass are allowed to vary with L (S-wave is quoted here), and a 
higher-mass a with large (;nNct couplings is also assumed. 

''Assumes rrip = m^^ . 

'^This value is cited in Ref. ll2l but is not actually used in the Paris 
model, which does not incorporate p exchange. 



